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Abstract 

A high precision numerical analysis of the static, spherically symmetric SU (2) magnetic 
monopole equations is carried out. Using multi-shooting and multi-domain spectral meth- 
ods, the mass of the monopole is obtained rather precisely as a function of /3 = Mh/Mw 
for a large /3-interval {Mh and Mw denote the mass of the Higgs and gauge field respec- 
tively). The numerical results necessitated the reexamination and subsequent correction of 
a previous asymptotic analysis of the monopole mass in the literature for /3 <C 1. 



1 Introduction 

It has been found some time ago by 't Hooft and Polyakov [1, 2] that spontaneously broken 
gauge theories (with a simple gauge group) admit classical solutions with quantized magnetic 
charge and finite energy. The self-energy (or mass) of these magnetic monopole solutions is well 
localized at any instant of time, hence they can be interpreted as 'particles'. Furthermore these 
magnetic monopoles are stabilized by a quantum number of topological origin, corresponding to 
their magnetic charge. Magnetic monopoles provide a natural explanation of the quantization 
of the electric charge, and they inevitably arise in grand unified theories. They have been very 
extensively studied, see the recent monograph of Manton and Sutcliffc [3] for a very detailed 
description of the zoo of topological solitons, including magnetic monopoles and a bibliography. 

An interesting problem for magnetic monopoles (or more generally for particle-like so- 
lutions of finite energy) is a quantitative description of their internal structure, their size and 
their mass in function of the parameters of the underlying theory. Since the equations describ- 
ing such objects are nonlinear, simple analytic descriptions are rarely available and in most 
cases one has to resort to numerical techniques (often combined with some analytic methods). 
The simplest prototype of the 't Hooft-Polyakov monopole arises in an SU(2) gauge theory cou- 
pled to a Higgs scalar in the adjoint representation. There arc two mass scales in this model, 
the mass of the vector mesons, Mw, and the mass of the Higgs particle. My. The monopole 
solution of 't Hooft and Polyakov is static, spherically symmetric and it depends on a single 
dimensionless parameter, (3 = Mh/Mw- An analytic solution is available in the special case 
/? = [4, 5], for non-zero values of (3 monopole solutions have been studied mostly numerically 
e.g. [1, 6, 7] and [3]. In the two limiting cases, /3 <C 1, and /3 2> 1, a combination of asymptotic 
and numerical techniques [8, 9] have permitted to derive asymptotic expansions for the mass 
of the monopole as a function of the parameter /3. 

The aim of this paper is to adapt and test numerical techniques based on combining 
standard Rungc-Kutta integration with multi-domain spectral methods, which latter have been 
successfully applied for 3 dimensional elliptic problems in Numerical Relativity [10]. Our results 
demonstrate that one can achieve very good numerical precision for a very wide range of the 
parameter (3, 10~^ < /5 < 2 x 10^, and we tabulate the values of the pertinent parameters of 
the solutions to 11 significant digits. By comparing our high precision numerical result with 
the asymptotic expansion of the monopole mass for /3 <C 1 presented in Ref. [9] 



we have identified a discrepancy for the coefficient C3, in that we have found C3 ~ 0.4429 from 
a numerical fit. This discrepancy has led us to reexamine the computation and we have found 
a computational error in [9] . Moreover we have found the following simple analytic expression 
for the coefficient C3: 

^ 2 24 72 

We have also determined the first three coefficients of the asymptotic expansion of the monopole 
mass in the limit /? — 00, confirming and even improving on the previous numerical results of 



We believe that the numerical techniques developed in this paper will prove useful to 
obtain good numerical precision for a much larger class of nonlinear equations, where simpler 
methods are not readily applicable. 




C3 = 0.7071 



[8]. 
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2 The Monopole Equations 

The action of an SU (2) Yang-Mills-Higgs theory with the Higgs field in the adjoint represen- 
tation is given as 

5ymh = j d'x (^-^F^^F'^'^^ + i(D^$'^)(L>'^$«) - |($»$« - v'^f^ , (1) 

where a = (1,2,3), e is the gauge couphng constant, 7 resp. v denotes the couphng of the 
scalar self-interaction, resp. the vacuum expectation value of the Higgs field. The gauge field 
strength and the covariant derivative are defined as 

C = ^M^" - + e^'^A^Al , L>^$« = 0^$" + e'^^M^^'^ . (2) 

The mass of the Higgs, resp. of the gauge field is given by Mh = v^/^, resp. Myy = v e. The 
'minimal' static, purely magnetic, spherically symmetric Ansatz for the gauge and Higgs fields 
can be written as 

^g(r) = , 4"(r) = 6,„fe J (Wir) - 1) , $"(r) = H{r)^ , (3) 

where x" are the cartesian coordinates and = x'^Xk- This Ansatz leads to the following 
reduced action: 

5™„ = _i= fw- + + + ^iH'- If- + W'hA . (4) 







e 

where the coordinate r and the function H have been rescaled as 

r 



ev 



, H^vH, and P = Mh/Mw (5) 



The equations of motion are obtained by varying the action (4), and they can be recast as a 
non-autonomous dynamical system of the form 



W'=U (6a) 

u' = (6b) 

H' = V, (6c) 



2 r 

where W, U, H and V are now considered as independent phase-space variables. 

Since the configuration is static, the self-energy, i.e. the mass of the monopole is given 
by = —Symh, and to simplify matters we shall use the rescaled mass, E, 

~ Attv 
E = MoE{P) , where Mq = . 

e 

Since we shall study in some detail the (3 dependence of the mass function of the monopole, 
it will prove useful to compute the derivative of E{P) w.r. to /?, and one finds: 



dEjp) _ dEjP) _ (5 
dp ~ d/3 



^J^^drrHH'-lf>0, (7) 



where the equations of motion (6) have been used. Equation (7) implies that the mass of the 
monopole increases monotonically with increasing /3. 



Let us recall that for the case of a massless Higgs field, i.e. /? = 0, the solutions of Eqs. 
(6) are known in closed form [4, 5] 



Wo{r) 



sinh(r) 



ifo(r) = coth(r) 



(8) 



which is the celebrated Bogomolnyi-Prasad-Sommerfield (BPS) monopole solution. In fact the 
BPS solution satisfies the following equations: 



0-"0 ) 



H' 
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(9) 



the spherically symmetric Bogomolnyi equations [4], which are much simpler than the original 
field equations (6). The (rescalcd) mass of the BPS monopole is just £'(0) = 1. 

In the limit /? — > oo the potential energy term in (1) forces the Higgs field to be frozen 
at its vacuum value almost everywhere {H{r) = 1 Vr > 0, but H{G) = 0) and the equations of 
motion (6) reduce to a massive Yang-Mills equation for the limiting solution (r) 



W" 



The mass of the monopole stays finite in this limit, and it is given by 



E{oo) 



dr 



2r2 



(10) 



(11) 



3 Numerical Techniques 

We briefly recall the local behaviour of the functions H and W ioi r ^ resp. r ^ oo, 
corresponding to the two explicit singular points of the dynamical system (6a-d). The boundary 
conditions at these points follow from the requirement of regularity of the solutions, which 

ensures the finiteness of the energy. 

One can show that there exists a two-parameter family of local solutions regular at r = 
given by the expansion: 

W{r) = l- br^ + 0{r'^) , H{r) = ar + 0{r^) , (12) 

(as long as /? is finite) with free parameters a, b. In fact from the finiteness of the energy (4) 
alone one can already deduce that 

H{r) 1 , W{r) ^ , for r ^ +oo . (13) 

In order to obtain globally regular, finite energy solutions of Eqs. (6a-d) by numerical inte- 
gration from the origin, one has to suppress the exponentially growing modes present in Eqs. 
(6a-d) by finc-timing the parameters a,b defined at r = 0. To exhibit the exponentially growing 
modes in the equations of motion for /? / 0, we introduce the following new variables 

W+ = W + U , W- = W-U , h+ = (3h + v , h- = /3h-v, (14) 

with h = r{H — 1) and v = h', which satisfy the following dynamical system 

W!^ = W+ + Fw{W, h, r) , h'^ = f3h+ + Fh{W, h, r) , 

(15) 

W'_ = -W- - Fw{W, h, r) , h'_ = -ph^- Fh{W, h, r) , 
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where 

Wh h WiW'^ - V) 

W,/.,r) = ^^(l + ^) + ^(3 + ^). (16b) 

Prom Eqs. (15) it is clear that W+ correspond to exponentially growing modes as r increases, 
while h-, W- are decreasing modes. Finiteness of the energy requires the absence of the growing 
modes for r oo, i.e. the functions, hj^, W+, should tend to zero for r — ^ oo. This implies 
that the asymptotic behaviour of iJ, W regular for r — > oo is given as 

W{r) = C,,e-'[l + 0{^)^ , (17a) 
H{r) = 1-Ch ^ [l + 0(^)^ , for /3 < 2 (17b) 

2(72 -2r / 1 \ -/3r / i \ 

Hir) = l-^^^[l + 0{-)yC,^(l + 0{-)y for/3>2, (17c) 

where Cu, Cvy &re free parameters determined only by global requirements like regularity at 
r = 0. Note that for (3 > 2 the asymptotic behaviour of H is dominated by that of W induced 
by the nonlinear terms (compare Eq. (17c) with Eq. (27) of Ref. [8]). For large values of (3 the 
two divergent modes, i.e. that of the Higgs field ~ e^'', and that of W (~ e''), are quite different 
that is we have a stiff system. 

Let us now sketch the numerical methods we have adapted. Since the locally regular 
solutions at the singular point r = admit a power series expansion (12) with two undetermined 
parameters a and 6, it is perfectly adequate to integrate the equations (6a-d) from the origin 
using standard Runge-Kutta methods. Actually, we have implemented an adaptive stepsize 
fifth-order Runge-Kutta algorithm [11]. 

At r = +00, the locally regular solutions do not admit, however, a power series expansion, 
therefore we have chosen to use integral equations, which are well suited to deal with such 
singular points. The solution of the system (15), satisfying the condition (13) is given by the 
following integrals 

W+{r)= j e''-''' Fw{W,h,r')dr' , (18a) 

J oo 

W-{r) = W-(rm)e''^~'' - [ e^'"'' FwiW, h, r') dr' , (18b) 
h+{r) = r e^('-''') Fh{W, h, r') dr' , (18c) 

J oo 

^_(r) = /i_(r„)e^(''--^) - T e^^^'"^) F^(T^, /i, /) , (18d) 

for some Tyji > 0. In order to solve these integral equations numerically we have compactificd the 
'outer' region, [rm, +oo), by introducing the variable t = 1/r. Then, the interval [0, tm = ^/rm] 
is discretized and the solution is obtained by iteration using 

W+{r) = 0, W-{r) = W-{rm)e''^-''' , 

(19) 

h+{r) = 0, h-{r) = /i_(r„,)e^('^--'') . 

as an initial configuration. Our general strategy to find globally regular solutions is to fine-tune 
a and h together with W-{rm) and h-{rm), in order to match the functions W, U, H and V 
at r = Tjn- 
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The above procedure works well and yields accurate results for 10~^ i$ /5 ^ 7. For 
values of (3 < 10~^ one has to increase more and more to ensure convergence of the it- 
eration of the integral equations (18a-d). This leads to a loss of accuracy, however, mainly 
due to the increasing length of the integration interval. To overcome this problem we have 
implemented a suitable version of the 'multi-shooting' method, by introducing some inter- 
mediate domains between the inner region (containing the origin) and the external one 
{[0,tm = l/rm])- We have integrated in each domain independently and then performed a 
global matching to obtain a regular solution by fine-tuning the set of 4 x (m -|- 1) free parame- 
ters, {a,b, W- {r„i) , /i- (r-m) , W{ri) , U{ri) , if (r.j) , F (r,) . | i = 0, . . . , m — 1}, where m denotes 
the number of domains. For small values of /?, the dominant divergent mode comes from the 
gauge sector ~ e*", and the choice of intermediate regions of roughly equal size ~ 7 has turned 
out to be adequate. 

For large values of /3, the situation changes radically, since the dominant divergent mode 
then comes from the Higgs sector and grows as '-^ e^''. Runge-Kutta algorithms fail to cope with 
such exponentially growing modes on intervals of size S> Therefore one should introduce 
m ~ /? domains for large values of /3, leading to a huge increase in CPU time and to potential 
numerical problems such as instabilities. To overcome this problem we have implemented a 
simple version of multi-domain spectral methods which are of global nature and are known 
for their high accuracy, see Ref. [10] for their use in Numerical Relativity. These techniques 
have turned out to be quite suitable to deal with such exponentially growing modes. In the 
procedure we have implemented, the values of W+, h+ and /i_ are iteratively computed 
starting from an initial guess. In each domain we parametrize the solutions by the values of 
W- and /i_ at and and /i+ at rj+i in order to control the growing modes during the 
iteration. We describe briefly some of the basic ideas of spectral methods in Appendix 1. 

The two cases, (3 = and /? = oo, require some special considerations. When /3 = the 
asymptotic behaviour of the BPS solution is dominated by the long range interaction of the 
massless Higgs field and the generic behaviour of H at infinity is 

Hir) = do + - + o(r-2) . (20) 
r 

where do, di are free parameters. The suitable combinations in the Higgs sector are then given 

by 

h+ = H-l + rV, h_=r'^V, (21) 

and the BPS solution (8) is obtained by imposing the boundary conditions = = at 
infinity. 

In the limit P — > oo, one has to integrate the massive Yang-Mills equation (10). Due to 
the singular behaviour of the Higgs field for P = oo, the function remains differentiable 
but not at the origin 

W^{r) = 1 + Hr) - b'oo + 0{r^ ln(r)) . (22) 

The behaviour of at infinity is not modified and one has to solve the integral equations 
(18a-b) with Fw{W^, h = 0,r). 

4 Numerical Results 

We present the numerically obtained vahies for the parameters a and b as well as for the energy, 
E, in Tables 1, 2 and 3. Figures 1 show H and W, while Figure 2 displays the mass function, 
E{P) for < /3 < 50. We remark that our numerical data give strong evidence that H resp. W 
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are monotonically increasing resp. decreasing functions of the radial variable, r, for all values 
of p. 

We expect that our numerical results for a, E, are accurate up to 11 significant digits. 
The parameter b is somewhat less accurate for large values of /3: its accuracy is around 9 
significant digits. The expectation for the quoted accuracy is based on the stability of the 
results w.r. to any change of the control parameters of the numerical integration. We shall 
use our high precision numerical results to extract the asymptotic behaviour of the monopole 
mass, E{I3) for /? — > and for /3 



oo. 





Figure 1: Plots of H and for /3 = (solid), 5 (dashed) and 50 (dotted). Note the difi^erent 
scales for H and for W. 




Figure 2: The monopole mass, E{P), plotted against (3. 

We study first the behaviour of the monopole mass for values of /9 <C 1. The analytically 
known value of the energy for /? = 0, -E'(O) = 1, provides a good test for the numerical results. 
In fact we have obtained numerically E(0) = 1.0 to an accuracy better than lO"^-*^. As it has 
been shown in Ref. [9] for /3 <C 1 the monopole mass, E{P), can be expanded as 

E{P) = Co + ci /3 + C2 /3Mn/3 + cg /^^ + o(/?2 ln/3) . (23) 

To determine the coefficients q from our results, we have computed E{(3) for 81 values of (3 
in the range 10"*^ < /3 < 5 x 10~^, and performed a standard least square fit, implementing a 
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X^-fitting routine as described in Ref. [11]. We have first minimized to obtain an estimate 
for the values of the coefficients Cj, and then determined the corresponding mean standard 
deviations, and cTc^ of E(P) respectively of the estimated values of the coefficients q. This 
way we have obtained 

Co = 1 ± 8 X 10"" , ci = 0.4999992 ± 1.2 x lO""^ , 

C2 = 0.496 ± 2 X 10"^ , C3 = 0.41 ± 1.5 x lO'^ , a^, ~ 3 x 10"" . (24) 

The values of the coefficients co, ci, C2 arc all in good agreement with the theoretical results 
of [9, 8]. The value, C3 = 0.41, however, deviates quite significantly from the one found in Ref. 
[9], where the value of C3 = 0.7071 has been quoted. In order to improve upon the precision of 
the value of C3 we have also made a fit on 

E{P) = 1 + 1/3 + 1/32 In /3 + 03/3^ + ^/33 ln/3, (25) 

i.e. fixing cq, ci, C2 at their theoretical values, and moreover we added an extra term, d (3^ In l3, 
to parametrize the contribution of the unknown higher order corrections. This way we have 
obtained 

C3 = 0.44297 ± 1.8 x 10"^ , d = 0.64 ± 5 x lO'^ , cj^ ~ 3 x 10"^^ , (26) 

a value fully consistent (within 2(Tc3) with the one found in Eq. (24). Without the "improvement 
term" , d(3^ln /3, the coefficient in question is given as C3 = 0.440890ib3.4x 10~^ . Being confident 
as to the correctness of the numerical results, we have been led to reconsider the analysis of 
Ref. [9]; this will be outlined in the next section. 

Next we consider the asymptotic behaviour of £^ for /3 3> 1. Motivated by the results of 
Ref. [8] we postulate the following asymptotic expansion for E: 

^ = e. + ibi + ^ + ^ + „(,-3). 

Following the same procedure as for /3 ^ 1, we have determined the numerical values of the 
coefficients e_j from a fit on 21 values in the range 1000 < /3 < 2000 : 

eo = 1.786658424 ± 1.7 x 10"^° , e_, = -2.228956 ± 7.2 x 10"^ , 
e_2 = 7.14±lxl0-2, e_3 = -32.7 ±4.5 (7^-3x10"". (28) 

These results provide some motivation to prove the validity of the expansion (27) by analytic 
methods. The value of £^(00) (corresponding to Cq) has been obtained independently from Eq. 
(11) and is given in Table 3. Note the good agreement between these two values. 

Some comments are in order on the obtained precision for the values of a, b, and on their 
respective /^-dependence. Since for (3^1 the Higgs field rises steeply in an interior region 
of size ~ 1//3 close to its vacuum expectation value (see also Fig. 1), H approaches more and 
more a step-function as /3 increases. Hence, in order to numerically integrate the equations 
near the origin, one has to choose an interior region of size ~ 1//3. Such a choice enables us 
to determine the parameter a with a good accuracy. Indeed, H is a linear function of r within 
some interval [0, C//3], therefore we expect the dominant term of a to be a linear function of 
/3. Indeed, taking 

a(/3):^ai/3 + ao + ^ + ^, (29) 

and fitting the coefficients we have found: 

ai = 0.357826169 ± 2 x 10"^ , Oo = 0.000677 ± 7.5 x 10"^ , 
a_i = 6.37 ± 10-2 a_2 = -440 ± 5 , (7^ ^ 8.2 x lO"*^ , (30) 
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for values of (3 ranging from 1000 to 2000. 

Related to the above choice for the interior region there is, however, some loss of precision 
on the value of b. This is due to the fact that W, as opposed to H, changes very little near 
the origin. We have found that the precision we get for b is only 10~^. This loss in precision 
for the parameter b does not significantly affect the precision obtained for the energy, which is 
an integrated quantity. Concerning the behaviour of b for large /3, we would like to point out 
that b'^ as defined in (22) does not correspond to the limit of 6 as /? — >■ oo. That is why the 
value b'^ = 0.4843164140 is absent from Table 3. Indeed for /3 — > oo the parameter b diverges 
(similarly to a), leading to the ln(r)/3 term in (22). However, we expect its divergence to 
be weaker than that of a. We found the behaviour of b to be well approximated by 

b{/3) ~ (0.33333270 ± 2 x 10"^) ln/3+ (0.065700 ± 2 x 10^^) + (6.82 ± 1.8 x 10"^) (3''^ , (31) 

with (Tft ~ 3 X 10~^. We note that the numerical value of the coefficient of the ln/3 in (31) 
is quite close to 1/3 which is precisely the coefficient of the ln(r) term in (22). It is very 
natural to conjecture that 6(/3) ~ (1/3) In/? as /3 — oo. 

Table 1: Numerical values of a, b and E for < P < 10~^. N denotes 
the number of domains of integration including the interior region and 
the exterior compactified region containing r = oo. 



f3 


a 


b 


E 


N 





0.3333333333 


0.1666666667 


1.0000000000 


2 


0.0001 


0.3333999248 


0.1666999621 


1.0000499584 


5 


0.0002 


0.3334663934 


0.1667331954 


1.0000998473 


5 


0.0003 


0.3335327549 


0.1667663745 


1.0001496747 


5 


0.0004 


0.3335990174 


0.1667995035 


1.0001994446 


5 


0.0005 


0.3336651866 


0.1668325852 


1.0002491600 


3 


0.0006 


0.3337312672 


0.1668656218 


1.0002988231 


3 


0.0007 


0.3337972625 


0.1668986152 


1.0003484357 


3 


0.0008 


0.3338631756 


0.1669315669 


1.0003979993 


3 


0.0009 


0.3339290091 


0.1669644781 


1.0004475152 


3 


0.0010 


0.3339947653 


0.1669973500 


1.0004969846 


3 



5 Asymptotic analysis of the monopole mass 

In this section we study the mass of the monopole in function of the mass ratio /3 = Mh/Mw 
for the two extreme cases /3 <C 1 (i.e. near the BPS limit) and for (3^1 (i.e. near the massive 
Yang- Mills limit). 

We present first a calculation of the monopole energy for /3 <C 1 to order This 
computation has been first performed in [9] based on matched asymptotic expansion techniques. 
The discrepancy between our numerical results presented in Section 4 and those predicted in 
[9] has led us to repeat the computation of [9]. 

It is natural to expand the monopole solution for /? ^ 1 around the BPS solution as 
= Hq + 5H , Wn = Wq + 5W and try to determine 5H , 5W , by linearizing the field 
equations (6a-d). The small perturbations 6H, 5W should satisfy re^iuiar boundary conditions 
at the origin. This approximation is valid, however, only in an interval < r < /3^"^ with 
some CKi < 1 ('near' region), but breaks down in the 'far' region, r S> /?~"^, where the Higgs 
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Table 2: Numerical values of a, b and £' for 10 ^ < /3 < 7. For this 
range of values of /3, only 2 domains are needed. 



/3 


a 


b 


E 


0.01 


0.3396968769 


0.1698451903 


1.0048108392 


0.05 


0.3624466559 


0.1811436175 


1.0220791858 


0.1 


0.3878538429 


0.1936144524 


1.0410559503 


0.5 


0.5533233202 


0.2697154617 


1.1493759746 


1 


0.7318137699 


0.3409021837 


1.2377010395 


2 


1.0683224404 


0.4490865847 


1.3510046819 


3 


1.4003124276 


0.5320943027 


1.4236384997 


4 


1.7338856977 


0.5997114280 


1.4749928035 


5 


2.0701446040 


0.6567291699 


1.5134382751 


6 


2.4090593765 


0.7059674601 


1.5433692766 


7 


2.7503251769 


0.7492501033 


1.5673591187 



field decays exponentially, H = Hp ~ 1 + Ae~^'^ /r + . . . with some constant A, as opposed to 
the slow 1/r falloff of the massless BPS solution. There is, however, an 'intermediate' region, 
p-a2 <r< , < a2 < CKi < 1 where both approximations hold, and one can match i?N 
with Hy- (The actual values of the exponents, ai, 02, do not (and should not) matter for the 
computation. In [9] the values ai = 1/2, 02 = 1/4 have been chosen for concreteness.) This 
matching is possible using the remaining free parameters in i^N , VFn resp. -ffp ) W-p after having 
imposed regularity at r = 0, resp. at r = 00. 

We rewrite the field equations (6a-d) introducing the rescaled variable x = /3r, which 
will prove useful in the far region and also for the /? » 1 limit: 

hV W{W^ - 1) 



W = W{ — ] +^ -, (32a) 



H = 2«(5,i,;.^_l,)_f. (32b) 

where / = df /dx. 

In the far region, a; 3> 1, the field equations (32) decouple since both fields fall off 
exponentially. H-p{x) then satisfies (32b) with W = Expanding H-p{x) around the free 
massive solution as 

/7p = l + ^p£:!+^2 /iMe-2. + ^3 ^Mg-3x+ ^ (33) 

XX X 

where A-p is the (small) expansion parameter, one obtains a set of linear differential equations 
for fn{x) which can be solved recursively. To determine C3 one only needs in fact fi{x), which 
satisfies: 

A-4A + 3/i = -, (34) 

together with the boundary condition fi{x) — > as x — ^ 00. The solution for fi{x) is then 
given as 

f^{x) = J (e^^Ei(-3x) - e^Ei(-ar)) , with Ei(-x) = / dt — . (35) 

Now it is not difficult to determine Hp{r) in the intermediate region, where x = Pr <^ 1, 
containing Ap as a parameter to be determined order by order in (3, 

Ay ( ^ , 31n3 , \ , A^p Mf„, \ „ , r,2^ 



^F(r) = 1 + ( 1 + 1 + -f- ( y 3 In r + r ) + i?F + o(/?') , (36) 
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Table 3: Numerical values of a, 6 and E ioi P > 10. 



/3 


a 


b 


E 


N 


10 


3.7849823090 


0.8542424691 


1.6173762372 


3 


50 


17.964403588 


1.371251333 


1.7447435157 


4 


100 


35.825856171 


1.60118232 


1.7650569743 


4 


500 


178.92488493 


2.13720411 


1.7822288542 


4 


1000 


357.83277188 


2.36823357 


1.7844365846 


5 


1100 


393.61488618 


2.40000240 


1.7846379843 


7 


1200 


429.39707896 


2.42900524 


1.7848059066 


7 


1300 


465.17933308 


2.45568539 


1.7849480580 


7 


1400 


500.96163610 


2.48038744 


1.7850699484 


7 


1500 


536.74397876 


2.50338457 


1.7851756213 


7 


1600 


572.52635402 


2.52489700 


1.7852681115 


7 


1700 


608.30875644 


2.54510487 


1.7853497407 


7 


1800 


644.09118172 


2.56415738 


1.7854223163 


7 


1900 


679.87362645 


2.58217954 


1.7854872652 


7 


2000 


715.65608785 


2.59927709 


1.7855457295 


7 


oc 






1.786(3581230 


2 



with 



Bp = At 



-1 + 



3Af 



ln3 , ^ 
— + ln/3 + 7E 



The above explicit form of Bp will not be needed for the calculation of the monopole mass, 
only the fact that Bp does not depend on r. 

In the near region it is somewhat more difficult to determine the solution, because the 
field equations are coupled there. The Green function in the BPS monopole background can 
be explicitly computed (using the Bogomolnyi equations) see e.g. Ref. [12]. This remarkable 
fact has been well exploited in [9] to find the solution H]<{, VFn. We omit here the somewhat 
lengthy computational details, since we basically agree with the results of Ref. [9], and give 
directly the final result for the expansion of i^N in the intermediate region: 



1 6^ 

- + ^(31nr 

r 2 



0^ 

3) + ^ 



31nA+ - 



+ An(A,/?) + o(/32 



(37) 



where r < A < is a kind of 'cutoff', which can be arbitrarily chosen within the intermediate 

region, A-!sf{A,f3) is the matching parameter. We remark that in contradistinction to Ap in 
the development of Hp, ^n(A, /?) turns out to depend explicitly on the chosen limit of the 
intermediate region, A. The value of Ap is now easily found by matching the r-dependent 
terms in Hp{r), Eq. (36) and in H-t^{r), Eq. (37), yielding 



Ap = -P - -{InS) /?2 + 0(/33) 



Bp = 13 



l + ^(31n(3/3) + 7E-l) 



+ o(/32). 



The value of An (A, f3) is then immediately found from the above results. Note the nonanalytic 
dependence of Hp(r) and Hj<f(r) on /? = \/]3^ and the appearance of terms oc Inr, ultimately 
responsible for correction terms of the type /31n/3 to the monopole mass. 

Now to compute the monopole mass for small values of /3 we evaluate first dE/d(3 (7) 
from which one easily obtains E{f3) integrating w.r. to /3. By splitting the integral in Eq. (7) 
at some matching point rm, such that /3~"2 < rm < A, 



+ EpiP,rn 



(38) 
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one can see that to first order in /? the two pieces rm), resp. E'p{/3,rm) are given as 

E^{f3,r^) = ^ rdrr\Hl-l)\ E^{P,r^) = '^ H drr\Hl-l)\ (39) 

with Hy taken to order in Eq. (39). For the details of the computation of dE/dj3 we refer 
to Appendix 2. From the result (64) found in Appendix 2 we obtain the development of the 
monopole mass for small [3: 

m = 1 + i/3 + In + 1 (^ln(37r) -^-^^p^ + o{p') . (40) 

Next we consider the /3 — 00 limit, already analyzed in detail in Ref. [8]. As shown in 
[8] the mass function can be expanded as 

The physical interpretation of this expansion is that the Higgs field does not contribute to 
the energy in the limit (3 00. Now, to evaluate dE/d{l/P) at /3 = 00 we take directly the 
/3 — 00 limit of Eqs. (32) which amounts to a 'blow up' of the origin to an infinite interval 
corresponding to the near region, where W = 1. Therefore the equation for Hoo is simply 

Hoo = 2H^ (1 + \{Hl - 1)) - ^ • (42) 

Note that H^o corresponds to the interior of a global monopole. In the outside ('far') region 
(i.e. r > 0), i?oo = 1 and Woo satisfies the massive Yang-Mills equation (10). 
The first correction to the mass is of order 1//3, and it is given by 
dE{P) i-oo ^2 

d{l/0) /3=oo" 

In order to compute this integral, one has to determine H^oix) by numerical integration of Eq. 
(42). The regular boundary conditions at the origin resp. at infinity are 



dx^-{Hl-lf. (43) 



H{x) = a'^x + 0{x-^), resp. H{x) = I - ^ - ^ - + ■ ■ ■ + B , (44) 

where and B arc a free parameters. Note that the behaviour of H{x) as ,t ^ 00 is dominated 
by an (asymptotic) series in 1/x with all coefficients determined, the free parameter, iJ, is in 
the sub-dominant term, /x. The numerical value of the coefficient (43) is found to be 

dE{l3) 

WM 

which is in very good agreement with the value obtained from our fit for e_i. It also agrees 
with the result found in Ref. [8]. The value of = 0.3578262475 agrees also quite well with 
ai in Eq. (30). This is due to the fact that H ~ o'^ x = a'^(5r, therefore a{f3) ~ a'^ (5. 



= -2.2289514397 , (45) 

/3=oo 



6 Conclusion 



We have illustrated the usefulness of multi-domain spectral methods combined with traditional 

Runge-Kutta integration to obtain high precision numerical results for the venerable 't Hooft- 
Polyakov monopole. As an application we have extracted the asymptotic behaviour of the 
monopole mass for small and large values of Our numerical results for small values of /3 are 
in disagreement with a previous result [9]. Repeating the asymptotic analysis of Ref. [9] we 
have found a simple analytic expression for the monopole mass, which is in full agreement with 
our numerical results. For large values of /3 our numerical results agree with those of Ref. [8]. 
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Appendix 1 



We present here some basics of spectral methods used in our approach. The fundamental idea 
underlying all spectral methods is to approximate a function / as 

N-l 

c,P,(x), (46) 

j=0 

where the Pj{x) form an orthonormal basis of polynomials of degree < (N — 1) and where 
f{xk) = P^{xk) for a certain set of N values {xq, ...,xn-i), called collocation points. We have 
chosen Chcbyshev polynomials, Tj{x), for the basis, Pj{x). Since they arc defined in the interval 
[—1,1] by the relation Tj(cos(,T)) = cos(jx), we shall consider the spectral decomposition of 
functions defined in this same interval. Each Chebyshev polynomial, Tj, has j zeros in the 
interval [—1,1]. A natural choice of the N collocation points where the interpolating polynomial 
coincides with the values of the function /, correspond to the N zeros of Tn{x) which are 

X. = cos , ^ = 0,l...,iV-l. (47) 

The Chebyshev polynomials satisfy the following discrete orthogonality relation 

N-i ( i ^ j 

^Ti{xk)Tj{xk) = { N/2 i = jy^O (48) 




fc=o { N i = j = 

Defining the coefficients cj , j = 0,1, . . . ,N 

2 ^"^ 

. |M/(=^.)cos(^^), m 

k=o ^ ^ 

the function / matches its polynomial approximation on the collocation points according to 

N-\ 

f{^k) = ^ + Y.cjTjixk). (50) 
j= 



It is interesting to note that the expressions (49) and (50) can be computed numerically using 
FFT algorithm if is a power of 2. 

For the integration of a system of first order differential equations by spectral methods 
one needs the primitive. Now, let F be a primitive of / and Cj , j = 0, 1, . . . , A — 1 the 
coefficients of the spectral decomposition of F. It turns out that a very simple relation relates 
the coefficients Cj to the Cj 

Cj = j > , with C7V = . (51) 

^3 

As a primitive of /, the function F has a single undetermined coefficient Cq. This coefficient 
has to be fixed by the value of F at some point in the interval. Usually, the value of Co is fixed 
at one of the boundaries of the interval [—1, 1]. Indeed, we have 

N-l N-l 

F{-l) = ^ + Y.Cun{-l) = ^ + Y.Ck{-lf, (52a) 

fc=l k=l 

„ N-l „ N-l 

F{i)=^ + Y.Ckn{i) = f + Y.Ck. (52b) 

k=l k=l 
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Appendix 2 



In this Appendix we present the computation of the pertinent integrals (39) in detail. 

We evaluate first the contribution to dE{P)/df3 from the far region, £'p(/3, rm). Using 
Eq. (33) to order for iJp one finds that it is given by 



1 f°° 

^F(A^m) = 4^/ dxx\HF{xf-l) 



j2 poo 



1 + Af + 2AFe-yi{x] 

X 



/32 



1 . / ln3 1 
-^F Inxm - a^m + 2 - I 7E + — + 2 



+ o(Xin) , (53) 



where .Tjn = /3ri„ (rm denotes the matching point). Using now Ap = — /3 — /?^(31n3)/4, Eq. 
(53) can be written as : 



E'p{P,Xm) = /^Inxm - + ^ + /5 Q +7E + ln3 ) + o{xm) ■ 



(54) 



To compute the monopole mass to order in the near region, < r < A, i^N should be 
identified with the BPS solution, 



1 



H^ir) = HqM = coth r , 

r 

and the contribution to dE{f5)/df3 from the near region is determined by 

rdrr\Hl-lf = 
Jo 



3 13 , , r(5sinh(r) — r cosh(r)) 

4r + coth(r) + ^ ^ , — 

r 3 ^ ^ 3 sinh'^(r) 

-^/i(rm)-8l2(rn,), 



r=0 



with 



/•r 

Jo 



dr- 



sinh (r) 



^ ^ ^ dr /I 1 1 \ 



The integral /i(rin) is not difficult to compute: 



TT 



/i(rm) = /i(oo) + o(l/rm) = Y + o(V''ii 



In order to compute the last integral, hirm), we consider the following identity 
hirui) = lim / dxF{p,x) - / dxF{p,x) } , 

P-^O I ^0 J2rrn J 

1 



with 



F{p,x) 



X \2 X — 1 
Now, the first integral of the r.h.s. of (58) to order 0(plnp) is given by 

/ dxF{p,x) = lnr(p) + -(Inp- In 27r) +p{l - In p) --(Inp + ln 27r) , 
Jo ^ 2 

while the second integral in Eq. (58) yields to order 



r°° 1 

/ dxF(p,x) = ~-Eii-2prm) + 

J2rm ^ 



-i\ P- 



.0 1 



(55) 



(56) 



(57) 

(58) 
(59) 

(60) 



(7£; + ln(2prn,) + 0(pr„i)). (61) 
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Therefore ^ 

/2(rm) = 2 (Inr-m + 7E - Ihtt) + o{l/r^) . (62) 

Putting now the above formulae together leads to the result (up to o{x-ai)): 

EN{P,r^) = ^ (4rm - 41nr^ - C) = - /3ln^ - /3 (^-fE + ^ + ^ - Invr^ . (63) 

Combining the contributions from the far (54) and from the near (63) regions, one finally 
obtains 

^.^,,,,,(^,K3.)-i-g)..<.(«. («) 

Note: Equation (6.3) of Ref. [9] corresponding to our Eq. (63) contains an important error. In 
[9] the value of the constant, C, is given (only numerically) as C = 0.9667, whereas the above 

analytical result yields C = 3.1599 We have also verified the correctness of C = 3.1599 by 

a direct numerical calculation. 
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